Theory of Magnetic Order in Fei +y Tei- x Se : 



Chen Fang, 1 B. Andrei Bernevig, 2 and Jiangping Hu 1 

1 Department of Physics, Purdue University, West Lafayette, Indiana 47907, USA 
2 Princeton Center for Theoretical Science, Princeton University, Princeton, NJ 08544 

(Dated: November 8, 2008) 

We develop a local spin model to explain the rich magnetic structures in the iron-based supercon- 
ductors Fei+ y Tei- x Se x . We show that our model exhibits both commensurate antiferromagnetic 
and incommensurate magnetic order along the crystal a-axis. The transition from the commensurate 
to the incommensurate phase is induced when the concentration of excess Fe atoms is larger than 
a critical value. Experimentally measurable spin-wave features are calculated, and the mean-field 
phase diagram of the model is obtained. Our model also suggests the existence of a large quantum 
critical region due to strong spin frustration upon increasing Se concentration. 

PACS numbers: 



Superconductivity with critical temperature T c = 1AK 
was recently reported in the iron-selenide-telluride com- 
pound Fei+ y Tei- x Se a r'££'& with < x < 1. This dis- 
covery not only adds a new class of iron-based supercon- 
ductors to the multitude of already existing ones, but 
also provides a fresh angle to investigate the fundamen- 
tal physics of the Fe — As based superconductors^^. 
Fe\+ y Te\- x Se x , similar to the Fe — As based materials, 
has a PbO structure of square planar sheets of tetrahe- 
drally coordinated Fe atoms. The electronic band struc- 
ture of Fei+ y Tei- x Se x calculated by LDA is very close 
to that of the Fe — As based superconductors 9 -. It ex- 
hibits electron pockets at the zone corner and hole pock- 
ets at the zone center. Similar to the Fe — As based 
superconductors, it is believed that the magnetism in 
Fe\+ y Te\- x Se x plays an important role in forming elec- 
tron Cooper pairs. 

The parent compounds of the Fe — As based supercon- 
ductors exhibit stripe-type commensurate antiferromag- 
netic spin order 10 . However, the origin of the spin order 
has been theoretically controversial, with two very dif- 
ferent mechanisms leading to the same physical answer. 
One theory is based on Fermi surface nesting between 
the electron and hole pockets at the zone corner and cen- 
ter respectively^. This weak-coupling approach leads to 
a commensurate spin density wave (SDW) state at the 
nesting wavevector, as observed in the experiments on 
Fe — As based parent materials. The competing view 
is that, due to the geometry of As- mediated hopping, 
antiferromagnetic exchange exists not only between the 
nearest neighbor (NN) Fe sites, but also between next 
nearest neighbor (NNN) site o 12 i 13 i 14 i 15 i 16 . Moreover, the 
NNN coupling strength J2 is stronger than the half of 
the NN coupling strength Ji/2. The resulting J\ — J 2 
model produces magnetic physics consistent with the ex- 
perimental results. Although they lead to the same over- 
all prediction, the two theories rely on different mecha- 
nisms to describe the parent state's magnetism, and a 
conclusive test of either of them is needed. 

We believe that the recent neutron scattering data 
in Fei+ y Tei- x Se x 17 i 18 sheds new light on the origin of 
magnetism in Fe-based materials. A commensurate spin 



order along the a axis (see FigpQ) in FeTe was reported. 
This gives way to an incommensurate spin order along 
the same axis upon the introduction of excess Fe atoms. 
The commensurate spin order in FeTe is different from 
the one in the Fe — As parent compounds: the two or- 
dered wavevectors are rotated by 45 degrees with respect 
to each other. This experimental fact places a clear chal- 
lenge to theories based on Fermi surface nesting. The 
Fermi surfaces of Fei+ y Tei- x Se x are predicted, from 
LDA studies, to be similar to those in the Fe — As based 
materials: the electron and hole pockets are separated 
by a 2D nesting vector at (71", 0). Based on first-principle 
calculations, it was then predicted that FeTe should sup- 
port an identical spin ground state to that observed in 
Fe — As materials 9 . With the experiment falsifying this 
prediction, the spin order in FeTe cannot be, at least 
trivially, understood by a Fermi surface nesting mecha- 
nisms. 

In this letter, we show that the magnetic physics in 
the FeTe parent compound can be understood from 
the usual magnetic exchange nearest and next nearest 
neighbor J\ — J2 model used for the Fe — As based 
materials, but with a natural parameter extension that 
takes into account the monoclinic lattice distortion ob- 
served in these compounds. The lattice distortion in 
Fei+ y Tei- x Se x is different from the one in Fe — As 
based materials. The two lattice distortion directions 
form a 45-degree angle, just like the magnetic wavevec- 
tors in the magnetic ordering states of these two systems. 
Our extended J\ — J2 model can explain both the com- 
mensurate and the incommensurate spin order phases 
along the a-axis which have been measured in neutron 
scattering experiments. The commensurate to incom- 
mensurate phase transition takes place at a critical con- 
centration of excess Fe atom in Fe\+ y Tei- x Se x . Above 
the critical concentration, the incommensurate wavevec- 
tor is proportional to the square root of the concentration 
difference of excess Fe atoms. 

We start with the J\ — J2 model on the tetragonal lat- 
tice. Due to their proximity in temperature, we strongly 
believe that the lattice and magnetic transitions in de- 
based materials are physically related. Considering the 
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coupling between the lattice and magnetism, it is phys- 
ical to assume that a particular lattice distortion favors 
changes in the values of J\ and J 2 as follows: in the 
Fe — As compounds, J\ should be slightly more sensi- 
tive to changes in the angle of the Fe — As — Fe bond 
(As is out of plane) than J2. This is because the angle 
between two nearest Fe atoms is around 72 degrees and 
hence closer to 90 degrees than the one between two next 
nearest Fe atoms which is around 112 degrees. However, 
in Fe\+ y Te\- x Se x , J2 should be significantly more sen- 
sitive than Ji to changes in the angle between two near- 
est Fe atoms; the angle between two nearest Fe atoms, 
which influences Ji, is around 66 degrees, whereas the 
angle between the two next nearest Fe atoms is around 
96 degrees and hence much closer to 90 degrees. 

In the monoclinic lattice distorted phase, the extended 
magnetic Hamiltonian can be written as an in-plane near- 
est and next nearest neighbor Heisenberg model supple- 
mented by an out-of plane small antiferromagnetic cou- 
pling as well as an next next nearest neighbor term: 
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where n is the layer index. The Ji/s, defined in Fig. [I], 
are the magnetic exchange coupling parameters between 
irons in the a— b plane, and their values depend on the lat- 
tice distortion direction. If small, the added J3 next next 
nearest neighbor coupling, suggested by first-principle 
calculations 19 , influences the phase diagram only quan- 
titatively. We take J2 a > ^26, Jia > Jib and study 
the part of the phase diagram of the model for which 
(^2a,^2b) > and J\ a > 0. These values are naturally 
expected in Fei+ y Tei- x Se x , as shown in Fig.pQ. 

As the exchange along the c-axis J z > is not frus- 
trated, we only focus on in-plane magnetism. All the 
possible ground states are presumed to have a k z = tt. 
The classical ground state of the Hamiltonian can be 
obtained by comparing the energy of the following six 
states:(l) (?r, tt) antiferromagnetic (AFM) phase (AFM1 
phase) with E 1 = -J\ a ~ Jib + J2a + ^26 + 2 J 3 ; (2 ) (0, tt) 
AFM state (AFM2 phase) with E 2 = - J 2a ~ J 2b + 2 J 3 ; 
(3) Commensurate AFM along the a-axis(AFM3) with 
E3 = -Jia-J2a + Jib + J2b-2Js (see Fig©); (4) Incom- 
mensurate AFM along the 6-axis(ICB phase) with energy 
E A = Min(-J 2a +(Jia-Jib)cos((f)) + (J2 b +2J 3 )cos(2(l))) : 
to be minimized over 0; (5) Incommensurate AFM 
along the a- axis (ICA phase) with E§ = Min(J2b + 

JlaCOs((/)i) + JlbCOs((j)2) + (J2a + 2J 3 )cOs(0i + 2 )), to 

be minimized with respect to 0i and 02 ; (6) Incom- 
mensurate AFM along the ?/-axis (ICY phase), with 
E 6 = Mm(-^ + (J 2a + J 2b - Ji±±Ji±) C0S ^ 3 ) + 

Jscos(2(j)s) + J3), to be minimized over 03. The three 
incommensurate phases are depicted in Fig. [2]. 

The phase diagram of above model in the — j^- 
plane is plotted in Fig. [2]. For FeTe, we are interested 
in the AFM3 and ICA phases, which have been experi- 
mentally observed. The AFM3 phase exists when the fol- 
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FIG. 1: (a) Fe (noted as Fe 1 ) square lattice showing commen- 
surate AFM3 order along the a-axis. The magnetic exchange 
parameters Ji a , Jib, J2a and J26 are specified, (b) The mon- 
oclinic lattice distortion observed in neutron scatterin g 17 ^ 18 , 
(c) The magnetic exchange coupling between Fe 1 and an ex- 
cess Iron atom in the middle of the plaquette indicated as Fe 2 . 



lowing two conditions are satisfied: J\ a > J16 + 4J26 — 8J3 



and 



< T J ]_ a Tt where JL = J 2a + 2J 3 . The tmnSi- 
tion line between AFM3 and ICA phases is determined 
by equality in the latter condition. In the ICA phase, the 
spin angle difference of two next nearest neighbours Fe 
atoms 0i and 02, indicated in Fig J2] are given by 
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The incommensurate wavevector, 5 along the a-axis cor- 
responding to the definition irtf is S = ^ 1 2 t ^ >2 . If the next 
next nearest neighbor coupling J3 is large, its effect on 
the phase diagram is qualitative. It can be analytically 
shown that if J3 > J26/2, then the ICB phase disappears; 
while if J3 > J2a/2, then the ICY phase vanishes. 

We now phenomenologically justify the parame- 
ters Jia-, Jib-, J2a-> J2b and their capacity to induce a 
commensurate-incommensurate magnetic order transi- 
tion upon doping with excess Fe atoms denoted as Fe 2 
(experimentally introduced in FeTe). The excess Fe 2 
atoms are located at the mirror-symmetric site of the Te 
atom with respect to the Fe layer. This suggests that 
the additional Fe 2 directly couples to the four nearest 
neighbor Fe as shown in Fig. [It] with coupling param- 
eter AJ. In the AFM3 phase, regardless of the sign of 
A J, the excess Fe 2 spins align along the b-axis through 
the AJ coupling with the collinear spins along the b- 
axis. The A J coupling along the a-axis is then frustrated. 
The effective result of A J on the magnetic model can be 
described by changing the effective parameters Ji a , Jib 
and J 2a as Ji a (y) = Jia + 2/ A J, Jib(y) = Jib + 2/AJ and 
J2a(y) = J 2a + 2/ 2 AJ, where y is the density of the ex- 
cess Fe atoms. The change of J 2a is a second order in 
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FIG. 2: The phase diagram of the Ji(Ji a , Jib) — Ji^Jia, J 2b) — 
J3 model: (1) AFM1 is the familiar (7r,7r) antiferromagnetic 
order in a square lattice; (2) AFM2 is the (0, 7r) collinear AFM 
order along the x-axis observed in the Fe—As based supercon- 
ductors; it is possible in the present model when J\ a — Jib, 
indicated as the green line in the figure; (3) AFM3 is the 
collinear commensurate AFM order along the a- axis observed 
in Fei+ y Tei- x Se x ;(4) ICB is an incommensurate phase with 
the incommensurate vector along the b-axis; (5) ICA is an in- 
commensurate phase with the incommensurate vector along 
the a-axis which has been observed in Fei+ y Tei- x Se x 17 ; (6) 
ICY is an incommensurate phase that has a wavevector of the 
form (7T,q). 



y, which can be ignored for low concentration of Fe 2 . 
According to the phase diagram in Fig. [2], the excess Fe 
atoms can cause a commensurate-incommensurate tran- 
sition if AJ is antiferromagnetic, i.e. positive. From 
the phase boundary between AFM3 to ICA we can then 
determine the critical concentration of Fe 2 that would 
cause the transition from AFM3 to ICA: 



V( J la ~ Jlb)(Jla + 4J* a - Ji 6 ) - ( J la + Ji 6 ) 
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Away from the critical concentration, the incommensu- 
rate angles as a function Sy = y — y c are given by 
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Therefore the incommensurate wavevector close to the 
transition is given by 
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To predict experimentally observable quantities, we 
perform a spinwave analysis on the AFM3 state. Sup- 
pose that the system is in the classical AFM3 ground 
state, in which the spins are aligned as in Fig.pQ. This 
is the commensurate order recently observed in the neu- 
tron scattering experiments. By introducing the stan- 
dard Holstein-Primakoff bosons we obtain the spin wave 
excitations, 
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where ei(fc) = 2( J ia + J 2a - Jib 
k y ), 62(h) = 2J 2a cos(k x + K) 



and f(k) = (e 
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lk y). The explicit analytical expres- 



sions for the spin-wave dispersion spectra are unreason- 
ably long and will not be given here. We plot the 3-D 
spin- wave dispersion in Fig. [3]. Regardless of the values 
of parameters, a common feature of the spin- wave disper- 
sion is an almost-dispersionless line along (fc±7r, k) where 
two branches become degenerate with an energy around 
4*V(Jia - Jib + 2 J 2a - 2 J 26 )( JTa - Jib ~ 2J 2b + 4J 3 ). 
By comparing the spin-wave dispersions for differ- 
ent J3, we can also determine the value of J3. In 
Fig. [3], we compare the dispersion along (fc, — fc) for 
J3 = 0.1,0.4,0.8 respectively. The major difference 
lies in the higher-branch, which turns from a convex to 
a concave-shaped line. This feature can be tested in 
inelastic neutron scattering to determine the value of J3 
explicitly. 

The spin moment correction in large S limit is: 



A5 = E / ( b hKk)dk 

i=l,2 JBZ 



In the vicinity of AFM3-ICA transition line, Am is about 
25%. It is worth noting that Am does not diverge at the 
phase boundary because the ICA (ICB) phases can be ob- 
tained from the AFM3 state by a continuous rotation of 
the magnetization starting from zero at the phase tran- 
sition line, unlike in other phase transitions where the 
magnetic order wavevector direction changes 15 . 

We now discuss the influence of the observed magnetic 
order in the AFM3 phase on the electronic properties 
of the material. DFT calculations show that FeSe and 
FeTe have a very similar Fermi surface structure to that 
in Fe — As based materials. In the AFM3 state, the 
meanfield Hamiltonian can be written as H m f = Ho + 
H M , where H = J2ka (k)e(k)^ a (k) is kinetic energy 
that describes the band structure and Hm is the mean- 
field energy of the spin ordering, 
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FIG. 3: 3D spinwave dispersion for the parameter set of 
J2a J2b 1, Jia 2, J ib 0.94 and J 3 0.4, 
which gives a point on the AFM3-ICA phase boundary on 
the phase diagram. Also note that along the BZ boundary 
(k zb 7r, k), the two sheets are degenerate and almost disper- 
sionless. Below: Different dispersion at different parameter 
regions with common parameters Jia — J2b = 1. All the 
three dispersions are drawn at AFM3-ICA phase boundary 
given by (from left to right) ( J la , J lb , J 3 ) = (1,0.72,0.8), 
(Jia,Jib,J 3 ) = (2,0.94,0.4) and (Ji a , Jib, J 3 ) = (5,0.9,0.1) 
respectively. 



by increasing 5e concentration, the magnetic frustration 
increases. The frustration can lead to a close competition 
between the AFM3 and AFM2 phases which can result 
in a strong quantum critical region controlled by a quan- 
tum critical point between the AFM3 and AFM2 phases 
or an existence of a spin liquid stated, which is an in- 
teresting open question to study in future. This physics 
will be critical to understanding unconventional trans- 
port properties in the materials 21 at high temperature or 




FIG. 4: The Fermi surface of the system in the magnetic 
ordered AFM3 order taking Ao = 0.3£i. The dotted red lines 
crop out the (folded) Brillion zone. 



with A(q) = A Dj=o( 1 + ^(- 1 ) J )^( < 7 - Qj) where q = 
(f , f = if) in the AFM3 state. The resulting 
Fermi surface is given in Fig. [3] and remains gapless even 
at a considerably large order parameter Ao = 0.3£i. The 
reason is simple: the modulation vector is (7r/2,7r/2). 
This vector, and its multiples, can only couple two elec- 
tron or two hole pockets, but cannot couple electron and 
hole pockets together. This is in sharp contrast to the 
1111 system, where the magnetic order induces a partial 
gap at both the electron and hole Fermi pockets. This 
distinctive feature from 1111 or 122 systems may be de- 
tected using ARPES. 

It is also important to discuss what happens upon re- 
placing Te atoms by Se atoms in the parent compound. 
Experimentally, it was shown that superconductivity de- 
velops in FeTei- x Se x upon increasing x. Based on our 
model, the AFM3 phase is generated from strong cou- 
pling of lattice and magnetic degrees of freedom. The 
AFM3 phase must coexist with monoclinic lattice distor- 
tion which is not favored in the pure a—FeSe. Therefore, 



upon increasing the Se concentration. 

In summary, we constructed a magnetic ex- 
change model to explain the rich magnetic order in 
Fei+ y Tei- x Se x . The model exhibits both commensu- 
rate antiferromagnetic and incommensurate magnetic or- 
der along the crystal a-axis and describes the transi- 
tion from the commensurate to the incommensurate spin 
order upon increasing the concentration of excess Fe 
atoms. Our model can be explicitly tested by experimen- 
tally measurable spin- wave features. Our model also sug- 
gests an existence of a large critical region due to strong 
spin frustration upon increasing Se concentration. 
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